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Abstract 

Using the idea of the penguin box expansion we find an analytic ex- 
pression for e'/e in the Standard Model as a function of mt, m s (m c ) 
and two non-perturbative parameters B^ 2 ^ and B^ 2 ^ . This formula 
includes next-to- leading QCD/QED short distance effects calculated 
recently by means of the operator product expansion and renormaliza- 
tion group techniques. We also derive an analytic expression for the 
upper bound on e'/e as a function of \V C \>\, \V n h/V C To\, Bk and other 
relevant parameters. Numerical examples of the bound are given. 
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This year |]J we have analyzed the CP violating ratio e'/e in the Standard 
Model including leading and next-to-leading logarithmic contributions to the 
Wilson coefficient functions of the relevant local operators. Another next-to- 
leading order analysis of e'/e can be found in 0. 

Imposing the constraints from the CP conserving K — ► ttti data on the 
hadronic matrix elements of these operators we have given numerical results 
for e' j e as a function of Ajjjg, m t and two non-perturbative parameters B^ 2 ^ 
and Bg which cannot be fixed by the CP conserving data at present. 
These two parameters are defined by 

(Q 6 (m c )> = Bf /2) (Q 6 (m c )>J vac) (Q 8 (m c )) 2 = Bf /2) (Q 8 (m c ))< vac) , (1) 

where 

3 

Qe = (S a dp) v _ A (Q(3 Qa) Y+A Q 8 = -Z (S a ^) V -A C 1 ^)v+A ( 2 ) 

q Z q 

are the dominant QCD and electroweak penguin operators, respectively. The 
subscripts on the hadronic matrix elements denote the isospin of the final 
7T7r-state. The label "vac" stands for the vacuum insertion estimate of the 
hadronic matrix elements in question for which B^ 2 ^ = B^ 2 ^ = 1. The 
same result is found in the large iV limit || ||. Also lattice calculations give 
similar results B ( ^ /2) = 1.0 ± 0.2 g g] and Bf /2) = 1.0 ± 0.2 @, |, g |]. 
We have demonstrated in [|IJ that in QCD the parameters B^ 2 ^ and B^ 2 ^ 
depend only very weakly on the renormalization scale fi when fi > 1 GeV is 
considered. The \i dependence for the matrix elements in ([]]) is then given to 
an excellent accuracy by l/m%(/J<) with m s (/x) denoting the running strange 
quark mass. 

In the present letter we would like to cast the numerical results of [[|] into 
an analytic formula which exhibits the m t -dependence of e'/e together with 
the dependence on m s , B^ 2 ^ and B^ 2 \ Such an analytic formula should be 
useful to those phenomenologists and experimentalists who are not interested 
in getting involved with the technicalities of ref. M. Combining the analytic 
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formula obtained here with the analytic lower bound on m t derived recently 
by one of us we will be able to find an analytic upper bound on e'/e. 

In order to find an analytic expression for e'/e which exactly reproduces 
the results of [|TJ we use the idea of the penguin-box expansion (PBE) sug- 
gested and developed in [[10] . This method allows to express flavour-changing 
neutral current (FCNC) processes, in particular e'/e, as linear combinations 
of universal m t -dependent functions resulting from various penguin and box 
diagrams common to many FCNC processes. The m t -independent coef- 
ficients of these functions are fi and renormalization scheme independent. 
They depend in the case at hand only on A 
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The basic ideas of PBE have been discussed at length in ref. |10| and 
will not be repeated here. The simplest method for transforming the results 
of ref. |1[ into the analytic formula given below is presented at the end of 
section 3 of ref. flQ| . In this letter we give only the final formula for e'/e. 

The resulting analytic expression for e'/e is then given as follows 



e'/e = 10" 



ImAt 



1.7 x 10" 4 



F(x t ) 



(3) 



where 



F(x t ) =P + Px X(x t ) + P Y Y(x t ) + P z Z{x t ) + P E E(x t ) . (4) 



Next 



lm\ t = | Kb | |Kb| sin 5 = rj A 5 A 2 



(5) 



in the standard parameterization of the CKM matrix |TT[ and in the Wolfen- 
stein parameterization |T2|, respectively. Here 



and 



A = \V US \ = 0.22 



rj = Rb sin 5 



Rh — T 



Kb I = AX 2 



A 



K 
K 



ub 



cb 



(6) 
(7) 
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The basic m t - dependent functions are given by 



X(x t ) 
Y(x t ) 
Z(x t ) 



E(x t ) 



x t 



x t + 2 i 3x t — 6 
x t - 1 
x t -A 



(x t - l) 2 
St - 1 ' (x t - l) 2 



+ 



In 

In x t 



1 , 18xf - 163x, 3 + 259a; 2 - 108x 4 
'9 lnX ' + 144(x f - l) 3 + 



32xf - 38a; 3 - 15x 2 + 18z t 



72(x t 



lnxt 



--lnx 4 + 



x 2 (15- ISxt + Ax] 

6(1 -x t y 



In x t + 



Xt(18 — llxt — x 



12(1 -ar t 



(8) 



with Xt = rnf fM§f. In the range 100 GeV < m t < 300 GeV these functions 
can be approximated to better than 3% accuracy by the following expressions 





X(x t ) = 0.650 a# 
Z(x t ) = 0.175 x 



0.59 



0.78 



0.93 



Y(x t ) = 0.315x 
E(x t ) = 0.570 x~ - 51 . 



(9) 



The coefficients Pi contain the physics below the M\y-scale. They are given 



in terms of B^ 2 ^ 



B$' 2 \ 



m c ), B^ 2 ^ = 5g 0/Zj (m c ) and m s (m c ) as follows 



?(3/2), 



Pi 



(0) , 



150 Me V 

m s (m c ) 



{rfB^+rfB^). 



(10) 



The Pi are yU-independent and renormalization scheme independent. They 
depend however on Attf = A™.. In the table below we give the numerical 



MS 

values of rf\ rf^ and rf^ for different values of A 



jy[g at — ?77, c 
(0) „(6) 



1.4 GeV 

and the NDR renormalization scheme. The coefficients rf , r) DJ and rf do 
not depend on m s (m c ) as this dependence has been factored out. does, 
however, depend on the particular choice for the parameter B^ 2 ^ in the 
parameterization of the matrix element {Qq)q. The values given in the table 



correspond to the central value for B\ 



(1/2) 



6.7 ± 0.9 as determined from 
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the CP conserving data in ref. JT|. Variation of B^ 2 ^ in the full allowed 
range introduces an uncertainty of at most 18% in the rj column of the 
table. Since the parameters rf^ give only subdominant contributions to e'/ e 
keeping B^ 2 ^ and rf^ at their central values is a very good approximation. 

For different fi and renormalization schemes the numerical values in the 
table change without modifying the values of the Pj's as it should be. To 



this end also B^ 2 ^ and B^ IZ> have to be modified as they depend on the 
renormalization scheme and albeit weakly on \x. 

AS = 1 PBE coefficients for various A^jg. 



?(3/2) 





A m = 0.20 GeV 


= 0.25 GeV 


% 


jo) 

' i 


' i 


J«) 

' i 


J°) 

' i 


J6) 
' i 


' i 





-4.323 


9.388 


2.109 


-4.406 


10.660 


1.957 


X 


0.988 


0.016 





0.958 


0.019 





Y 


0.763 


0.078 





0.729 


0.086 





z 


0.235 


-0.014 


-10.072 


0.297 


-0.015 


-10.899 


E 


0.359 


-1.207 


0.411 


0.339 


-1.327 


0.462 






= 0.30 GeV 


Aj^g = 0.35 GeV 


i 


i 


J6) 
i 


J8) 
i 


JO) 

i 


J6) 

i 


J8) 
i 





-4.484 


12.038 


1.792 


-4.557 


13.573 


1.609 


X 


0.933 


0.022 





0.910 


0.026 





Y 


0.701 


0.094 





0.676 


0.103 





Z 


0.361 


-0.017 


-11.779 


0.427 


-0.019 


-12.737 


E 


0.319 


-1.449 


0.515 


0.299 


-1.579 


0.573 



The inspection of the table shows that the terms involving and 
dominate the ratio e'/e. The function Z(x t ) represents a gauge invariant 
combination of Z°- and 7-penguins ]10[ which increases rapidly with m t and 
due to < suppresses e'/e strongly for large m t 0, The term 
which also plays some role represents the contributions of (V — A) <g> (V — A) 
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QCD penguins. The positive contributions of X(xt) and Y(xt) represent to 
large extent the (V — A) ® (V — A) electroweak penguins. Finally E(x t ) 
describes the residual m t -dependence of the QCD-penguins which is only a 
very small correction to the full expression for e'/e. 

Combining the formula (H) with the analytic lower bound on m t from e^ 
derived recently in || we can find an analytic upper bound on e'/e. Indeed 
the lower bound of || corresponds to sin 5 = 1 at which ImAt is maximal. 
Since F(x t ) decreases with increasing x t we find 



10" 



R h A 5 A 2 
1.7 x lO" 4 



F((x t ) min ) 



(11) 



where according to 



2 A 2 \A 2 B K R h 



1.2 



1.316 



(12) 



with the renormalization group invariant Bk giving the size of the hadronic 
matrix element (K°\(sd)v-A{sd)v-A\K°) . Formula (|IID gives the maximal 
value for e'/e in the Standard Model consistent with the measured value 



of £k as a function of |Kt>|> |Kib/Kb|, B K , B, 



(1/2) R (3/2) 



m<= and A- 



MS' 



It 



should be noted that this formula does not involve m t directly as m t has 
been eliminated by means of the lower bound of 0. It is evident that the 
upper bound increases with increasing |Kb|, IK^/Kbl, 



B K , B^ /2 \ A m and 



with decreasing m s and B 



(3/2) 



In fig. HI we show (e'/e) max as a function of |Kb| for different values of 



(1/2) R (3/2). 
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,B ( 8 



300 MeV, mJm r 



150 MeV or 



B K , three choices of (B, 
equivalently m s (l GeV) = 175 MeV and (Kb/Kb) = 0.10. The range for |V^b 
has been chosen as in |§ in accordance with the analyses of [15, [17[ [18 
and the increased 5-meson life-time tb = 1.49 ± 0.03 ps |I9] 



B 



(3/2) 



1 corresponds to the leading 1/N result or central values 



of lattice calculations. B, 



(1/2) 



advocated in ref. [3DL Finally, B, 



2, B, 

.(1/2) 



(3/2) 

B, 



(3/2) 



1 is representative for the case 
2 corresponds effectively to 
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the first case with a smaller value of m s . Comparing fig. [I] with the most 
recent messages from NA31 and E731 collaborations pi] , |22f 



„ , ,. , f 23 ±6.5 • 10" 4 NA31 
Re (e' s) = { v ; , 13 

[ (7.4 ±6.0) • 10~ 4 E731 

we observe that whereas the results of E731 are fully compatible with the 
cases considered here, the NA31 data lie above the bounds of fig. [I] if \V c b\ < 
0.040. This is in particular the case for B$ 2 ' = B^ 2 ^ = 1 which we advo- 
cate. 

In fig. [I] we only show (s'/e) max > 0. We observe that for B K < 0.5 the 
upper bound is very low and for Bq 2 ' = B^ 2 ^ = 1 and \V c h\ < 0.039 the 
ratio e' I e becomes negative. This is simply the result of the high value of mt 



required by Ek (see eq. (|T2|)) when Bk and \V c b\ are low. For Bk = 0.7 ± 0.2 
obtained in the 1/N and lattice calculations, {e'/e) max is in the ball park 
of the E731 result. We note however that for | V cb | < 0.040, B K < 0.8 
and B^ /2) = Bf /2) = 1 one has {e'/e) max < 5 x 10~ 4 . Since \V uh /V ch \ is 



expected to be smaller than 0.10 p3L |2J we should indeed be prepared for 



(e'/e) max ~ few x 10~ 4 in the Standard Model if 1/N calculations and lattice 
results are the full story. This has been already emphasized in [|T], but 
the bounds in fig. 0(c) make this point even stronger. On the other hand 
we should keep in mind the strong dependence of the upper bound ( [TTD on 
B K , B^ 2 \ B^ 2 \ m s , | Kb | and | Kb/Kb I- Future improvements on these 
parameters may bring surprises. 

Finally it should be mentioned that below the upper bound (sin 5 < 1) the 
dependence of e'/e on Bk, |Kb| an d | Kb/Kbl for fixed m t is more complicated 
and different than in (|TT|) . e'/e increases with decreasing B^ and |V^b| at fixed 
|Kib/Kb| until the bound is reached. It increases (decreases) with decreasing 
|Kb/Kb| for tt/2 < 5 < it (0 < 5 < it/2) at fixed m t , B K and |Kb|- 

We hope that the analytic formula for e' j e and the analytic upper bound 
on this important ratio presented in eqs. (H) and (|TT|) , respectively, should 
facilitate the future phenomenological analyses once the values for |V^b|, 
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|Kib/Kb|, Bk, Bq , Bg , m s , Aj^g and m t have been improved and most 
importantly e'/e accurately measured. 
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Figure 1: {e'/e) max as a function of \V c b\ for |V^b/Kb| — 0.10 and various 
choices of B-&- The three (s'/e) max plots correspond to hadronic parameters 
(a) Bi 1/2 \m c ) = Bf /2 \m c ) = 2, (b) Bi 1/2 \m c ) = 2, B^ ,2 \m c ) = 1 and (c) 
Bi 1/2 \m c ) = Bf /2 \m c ) = 1, respectively. 
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